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We present an accurate determination of the refractive index of hydrogen, helium, H2O and neon 
up to 35 GPa at ambient temperature. The experimental method is based on the combination of 
two interferometric signals of the Fabry-Perot cavity containing the sample in the diamond anvil 
cell. The data are put in perspective through the variation of the molecular polarisability with 
density. Significant electronic changes are observed. Interesting possibilities of the method are also 
illustrated: the high precision to finely probe phase transition; the simultaneous measurement of 
the volume to obtain the equation of state of transparent media with around 1% accuracy. 

PACS numbers: 07.35. +k, 64.30. +t, 78.20.-e 



I. INTRODUCTION 

There exist considerable evidences that pressure can 
dramatically change the electronic properties of a sys- 
tem. On the atomic scale, the variation of the electronic 
kinetic and potential energies with volume tends to make 
a contraction of the atomic size and a delocalisation of 
valence electrons. Consequently, spectacular transfor- 
mations have been observed, such as rehybridization of 
atomic orbitals, like in carbon (sp 2 ) to (sp 3 ) across the 
graphite to diamond transition ; insulator-metal phase 
transitions, like in 0-^, Xe£ or Csl 3 -; the transformation 
of ice from an hydrogen bonded solid to an ionic solid 
in symmetric iceA. Most changes of the properties of 
a system under pressure are related to the microscopic 
electronic reorganization, and the determination of the 
refractive index over a large frequency domain is proba- 
bly the most direct method to probe the electronic struc- 
ture. The refractive index expresses the response of the 
distribution of charges to a perturbing electric field. In 
order to extract a detailed change of the charge distribu- 
tion of a system with pressure through the measurement 
of the refractive index, a large frequency domain has to 
be investigated. This is unfortunately very difficult at 
high pressure. However, already interesting electronic 
evolutions have been obtained from measurements of the 
refractive index performed in the visible, for instance in 
hydrogen with the red shift of an electronic transition 
that parallel the closure of the electronic gap (see refer- 
ences below). 

Also, an accurate determination of the refractive index 
under pressure has a clear metrological importance. In 
Brillouin spectroscopy, where the refractive index of the 
sample and of the pressure transmitting medium have 
to be known to extract the sound velocity from the fre- 
quency shifts. In the newly developed laser shock mea- 
surements on pre-compressed targets to extract the shock 
front velocity from the phase shift of the VISAR (veloc- 
ity interferometric system for any reflector) diagnostics^. 



Finally, to measure easily the thickness of the sample 
chamber in the diamond-anvil cell (DAC) with a single 
spectrometer, when the sample or the pressure transmit- 
ting medium is transparent. This can be directly used 
to measure the Equation of State (EoS) of the sample 
(see below), or help to control the experimental condi- 
tions, for instance by evaluating the strain of the studied 
sample. 

In the DAC, the tips of the diamond anvils that en- 
close the sample form a Fabry-Perot interferometer. This 
geometry have already been exploited to measure the re- 
fractive index of transparent samples at high pressure. 
Most of the experimental methods developed so far have 
been based on the measurement of fringes contrast (— , — , 
£, _Lil) . Other were based on a controlled sample thick- 
ness (— r^) ■ In one work«i£, the combination of the mea- 
surement of parallel white beam interference spectrum 
and convergent monochromatic beam inferference rings 
allowed the determination of both refractive index n and 
thickness t of the sample. We present here an improve- 
ment of this method. 

With this setup, we have also been able to measure 
the Equation of State of helium, by estimating the sam- 
ple chamber volume changes upon compression. Such a 
method has been previously described for H 2 We 
show here that with numerical photograph technology 
and the accurate determination of the refractive index, 
its accuracy becomes comparable to other techniques. 
Therefore, it becomes an interesting alternative to X-Ray 
methods (X-Ray diffraction and X-Ray absorption^) for 
lights elements, particularly light liquids, or to adiabatic 
compressibility measurement of Brillouin spectroscopj^. 

Three goals have mainly motivated the present accu- 
rate refractive index measurements on four systems, He, 
H 2 , Ne and H 2 0. First, H, He and H20 are the main 
constituents of planetary interiors and the determination 
of their properties up to planetary core conditions is mo- 
tivating a lot of effort. Second, He and Ne are the most 
commonly used quasi-hydrostatic pressure transmitting 
media in DAC. In both cases, the determination of the 
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refractive index should be useful to analyze some spec- 
troscopic measurements. Finally, they have simple elec- 
tronic structures with different electronic organisations 
and the present data set could be interesting to test the 
confidence of ab-initio calculations in predicting the op- 
tical properties. 

The experimental method and a discussion of the un- 
certainties is presented in section II. The refractive index 
measurements are presented in section III. The analysis 
of the data in terms of the changes in the polarizability 
with density is presented in section IV. Applications of 
the method to the determination of the equation of state 
and to the study of phase transitions are described in 
section V. 



II. EXPERIMENTAL PROCEDURE 
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FIG. 1: Experimental setup. The light emitted by the laser 1 
(He-Ne, 7 mW or Ar) is extended by LI (simple lens, / = 30 
mm) and focused on the sample by the objective L2 (/ = 50 
mm) ; interference rings propagate through L2, Ml, M2, L3 
(objective, / = 80 mm) and are recorded by the digital camera 
CI. On the right, the PRL setup (doubled YAG laser, Fl: 
illuminating optical fiber, C2: digital camera, S: spectrometer 
optical fiber) equipped with a x20 Mitutoyo objective, allows 
pressure measurement, imaging of the sample and white light 
interferences measurement. 



A. High pressure techniques 

In these experiments we used a membrane diamond 
anvil cell described elsewhere (— ) with a large optical 
aperture (± 30°) and equipped with 400 or 500 /im ta- 
ble diamonds. Samples have been loaded in a rhenium 
gasket at room temperature, by the use of high pressure 
loading techniques. In the case of H2O, ultra-pure wa- 
ter was inserted in the gasket hole using a syringe. As 
explained below, thick samples are needed for a better 
accuracy of the refractive index measurements. For this 
reason, double rhenium gaskets have been used in the 
highest pressure runs. These gaskets are prepared in the 
following manner: two rhenium gaskets are pre-indented, 
then polished on one side and glued together. With this 
method, the thickness of an helium sample was of the 
order of 50 [im at 200 kbar, 30 % more than with a sim- 
ple gasket. The pressure gauge was ruby, calibrated with 
the quasi-hydrostatic pressure scaled, that allows pres- 
sure measurement with an accuracy of the order of 0.5 
kbar. The absolute accuracy is the one of the ruby scale, 
i.e at maximum 2% over the pressure range investigated 
here. 



B. White light interferences 

The refractive index setup is modified from LeToullec 
et al.'s one (see Fig.^J. The white light interference pat- 
tern in transmission is recorded over the range 500nm to 
700nm with a double substractive spectrometer (Dilor). 
A small area, varying between 20 x 20 and 10 x 10 ^m 2 , 
was selected to perform the measurement. Under high 
pressures (P > 250 kbar), the contrast drastically de- 
creased out of a central flat region, because of the strain 
of the diamonds under pressure. The wavelength differ- 
ence between two intensity maxima : 

AA =2^? W 



where t is the thickness of the sample, n{\) its refractive 
index, was determined. In that wavelength range, dis- 
persion could be observed only in H2O samples; in that 
linear dependency of 1/n on A was assumed in 
order to fit Eq. ^ 

C. Monochromatic interferences 

The setup presented in Fig. allows the measurement 
of the reflective Fabry-Perot rings produced by the cavity 
between the diamonds, when illuminated by a convergent 
monochromatic beam. We used an Argon laser (A = 
488.0 nm) for helium samples 1 to 3, and a He-Ne laser 
(A = 632.8 nm) for H2O, H2, Ne and He samples. 

The interference order k at angle i is defined by: 

k = \lt{n 2 -s\n 2 i) 1/2 - (2) 
A 

if k is an integer, a minimum in intensity is observed. 

At the center of the pattern, the interference order 
ko = 2n(X)t/\, can be deduced from white light in- 
terferences measurement. The angular position i p of 
the p th intensity minimum has an interference order 
k p = [ko] — p + 1. Refractive index and thickness can be 
deduced from two fringes (p and p') angular positions: 

r, sin Zn /c„ sin %w 

- 2 = (3) 

V K p 

and 

e = L. . >' \ (4) 

4 sin i p — sin i p > 

Interference rings were projected and recorded by a CCD 
camera (CI, see Fig.^l. The radius R p of the rings on CI 
is thus measured. Because of optical aberrations, mainly 
caused by the objective lenses L2, the relationship be- 
tween R p and i p was non trivial. To correct these abbera- 
tions, a reference spectrum was recorded. This spectrum 
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consisted in fringes formed by an empty indented and 
drilled gasket, placed in the DAC. The adjustment of the 
optical setup was performed as it is described herebelow. 
First, the cell was centered and perpendicularly oriented 
on the laser beam. Then Ml, CI, L3, L2, LI and Dl 
were positioned and the calibration spectrum recorded. 
All theses elements were mounted on (x, y, z, 9) linear 
and rotation stages. After calibration, all the elements, 
except the cell, remained fixed. In the first experiments, 
the cell was regularly removed to record the white light 
interferences pattern. The best position of the cell (in 
the x direction) was determined, within ± 10 /im, by 
optimization of the rings pattern stability with respect 
to Y and X translations. This allowed a reproducibil- 
ity of the index measurement of An = ±5.10 -4 . Each 
digital Fabry-Perot pattern was automatically processed 
by a program that evaluates thickness and index for all 
available pair of rings. This procedure allows the analy- 
sis of numerous spectra, and again, optimization of the 
x position of the sample. Tentative error bars have been 
estimated for each measurement, taking into account two 
factors: first, the uncertainty on rings radius R p ; second, 
the scatter of the refractive index given by different cou- 
ples of rings (Eq. EJ. 

The precision of refractive index measurements is lim- 
ited by the following effects: 

- the precision of k measurement by white light inter- 
ferences. The relative precision on fco estimate is better 
than 10 -3 (Afc < 0.2) for weakly dispersive media (He, 
Ne and H2) and of approximately 2.10 -3 for dispersive 
media (H 2 0). This measurement must also be carried out 
in the same sample region than the region illuminated by 
the laser beam for monochromatic interferences, in order 
to prevent any error caused by non-parallelism of the di- 
amonds. 

- the number of Fabry-Perot rings. They must be as 
numerous as possible to statistically increase the preci- 
sion of the measurement. In these experiments, a mini- 
mum of four rings has been observed. For this purpose, 
the angular aperture of the optical setup was increased to 
25°. The number of rings is also proportional to the sam- 
ple thickness; therefore, thick samples have been used. 

- the elastic strain (cupping) of the diamonds. This 
is the main source of uncertainty of our measurements, 
that limits the pressure range of this method. This cup- 
ping can be first observed by translating the cell in y or 
z directions and looking at the modification of the rings 
pattern. It can be noted that these measurements con- 
stitute an accurate mean of studying the strain of the 
anvils. At the beginning of the compression, inverse cup- 
ping, caused by the gasket strength, is observed. Upon 
pressure increase, the diamonds become flat and then 
cupping occurs. It has been observed approximately at 
the same thickness of a simple rhenium gasket (40 to 45 
/im for H2O samples), and for pressures varying between 
100 and 200 kbar. When pressure increases, it perturbs 
the rings pattern out of a central flat zone, up to a point 
that makes measurements impossible. To reduce this in- 



trinsic limitation, thick gaskets can be used because they 
increase the pressure at which cupping appears. This 
elastic strain also makes the focusing diagnostic less pre- 
cise. 

- the pressure gradients in the solid phase. The ruby 
chip can not be directly located where the sample was 
illuminated, because it perturbs the Fabry-Perot inter- 
ferometric rings. It was located at ~ 60 /im from the 
refractive index measurement point. In separate runs we 
have estimated the radial pressure gradient in ice VII, in 
the same sample geometry; this gradient is negligible up 
to 200 kbar and reaches 0.33 kbar//tm at 350 kbar. In 
ice, the pressure was corrected by this factor. Actually, 
this pressure gradient is also an intrinsic limitation of the 
method, because the extend of the source of interference 
rings is 100 /im, if i < 25°; at high pressure, rings with 
i < 15° are analyzed in order to reduce this error. 



D. EoS measurement 

An analysis of numerical images of the sample taken in 
transmitting light allows the measurement of the surface 
S of the sample, with a precision of 1 %. The volume Vs 
is then obtained by multiplying S with the thickness of 
the sample, obtained at the same time as its refractive in- 
dex. The uncertainty on volume determination is mainly 
caused by the uncertainty on S. No corrections taking 
into account the cupping of the diamond was necessary 
for the measurements presented in Sec. V, but these cor- 
rections are in principle possible. We made the assump- 
tion that the amount of matter trapped in the sample 
chamber remained constant upon pressure increase. In 
that case, the EoS V(P) can be simply deduced from 
Vs(P) by choosing a reference point Vq(Pq). For a given 
run, the molar volume of the sample can be expressed as: 



V mol {P) = V(P) 



VmoK-Prcf) 
V(P Ie{ ) ' 



V{P) being the volume of the sample at the pressure 
P and V(P xc f) the volume of the sample at the refer- 
ence pressure P rc f . We noticed that the precision of the 
volume measurement increased when the thickness of the 
sample decreases. A similar observation has already been 
made by Evans and Silvera, 1998, — . This is probably 
because the geometrical imperfections of the gasket hole, 
extensively described 8 , are smaller for thin samples. 



III. REFRACTIVE INDEX DATA 
A. Helium 

The refractive index data obtained in 4 runs up to 210 
kbar are shown in Fig. During the same experimental 
run, the orders of interference of the rings of the fringe 
pattern have been followed between each pressure steps. 
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FIG. 2: Upper: Refractive index of fluid and solid He as 
a function of pressure at ambient temperature. The first 3 
runs have been performed at A = 488 nm, and the last one at 
A = 632.8 nm. The parameters of the fitted functions (dotted 
lines) are presented in Table U Lower: Difference, in %, 
between experimental data and the fitted function. 



FIG. 3: Same plot as Fig. [21 for neon. Experiments have been 
performed at A = 633 nm. 



The good agreement between the data of the different 
runs data indicates that the uncertainty associated to the 
white light determination of the order of interference is 
negligible. For all data shown in Fig. [21 the error bars are 
smaller than the size of the symbols. The difference be- 
tween the fitted index law and LeToullec measurements 
is smaller than 2.10~ 3 , except for the highest pressure 
points where LeToullec measurements were less precise. 



B. Neon 

One run has been performed with a neon sample, up 
to 27 GPa. The evolution of the refractive index of neon 
with pressure is plotted in Fig. The refractive index 
of Ne had been determined previously to a maximum 
pressure of 1 GPa from capacitance measurements 
The present data are in very good agreement with this 
previous determination, in the pressure range of overlap. 
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C. Hydrogen 
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The evolution of the refractive index of hydrogen with 
pressure is plotted in Fig. 31 

There is a systematic difference between our data and 
previous determinations, either the one by brillouin scat- 



FIG. 4: Refractive index of H2 at A = 632.8 nm. Inset: 
enlarged region in the fluid phase. 
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FIG. 5: Measured refractive index of water, ice VI and ice 
VII, at A = 632.8 nm, with the corresponding fitted functions 
(dotted lines). 



tering— (1.5 %) or the one by white light interference— 
(around 1 %, which is within their estimated error bars). 
The matching of refractive index of ruby and hydrogen 
has been observed (annulation of the contrast) at 35 ±10 
GPa. This provides an anchorage point to the refrac- 
tive index curve. At this pressure, the refractive index 
of ruby is expected to be 1.77—; this value is in good 
agreement with the present measurements and that gives 
a good confidence in our highest pressure runs. 



Compound 



Refractive index 



P range 
(kbar) 



n = 0.8034 + 0.20256(1 + p) 12763 
n = -0.1033 + (1 + P) 052 
n = 0.67 + 0.33(1 + 4.3 x p) - 076 
n = 0.9860 + 0.08578 p- 01953 
n = 0.949 + 0.06829(1 + 11.8 P) - 2853 
n = 0.95 + 0.1144 p - 3368 
Fluid H 2 n = 0.900 + 0.4323(1 + 0.1769 P) 0164 
Ice VI n = 1.425 + 0.00255 P 

Ice VII n = 1.175 + 0.2615(1 + 0.101 P)°- 222 



Fluid He 
Solid He 
Fluid Ne 
Solid Ne 
Fluid H 2 
Solid H 2 



0.8-115 
117-202 

7-47 
50-270 

3-49 
52-350 
0.5-12 
12.7-22 
30-354 



TABLE I: Fitted forms of refractive index vs. P (in kbar) 



D. Water and ice 

For H2O, 5 experimental runs have been made, up 
to 35 GPa. Due to the dispersion of the refractive in- 
dex of H2O, of the order of 1 % over the 500 - 800 nm 
range, the determination of the order of interference from 
white light measurements is less precise and consequently 



small deviations between different experimental runs can 
be detected in the fluid phase. However, our measure- 
ments agree very well with precise interferometric mea- 
surements in the fluid phase — . At higher pressures, 
our data are within the scatter of a previous Brillouin 
scattering determination — . At pressures higher than 20 
GPa, the strength of ice causes pressure gradients in the 
sample, that can perturb the pressure estimation. To 
correct it, the pressure measured by the ruby chip Pr 
was then corrected by a factor AP = d x dP/dr, d be- 
ing the distance between the ruby chip and the center of 
the sample, and the last factor has been estimated in a 
separate run by mapping the sample chamber with ruby 
(in kbar//im): dP/dr = (Pr x 0.425) /i. The estimated 
error bars plotted in Fig. increase drastically as pres- 
sure increases, showing the decrease of the quality of the 
fringe pattern associated to the strain in the diamonds. 
Refractive index matching between ice VII and ruby has 
been observed at P = 40.5 ± 10 GPa. At this pressure, 
the refractive index of ruby is approximately 1.765—, 
1.2 % lower than the extrapolation of our measurements. 
That is within the uncertainty of our highest pressure 
measurements. 

The functions obtained by a fitting procedure to our 
experimental data are gathered in table [I] The use of 
power laws is justified, to a certain extent, by the fact 
that the Lorentz-Lorenz factor defined in Eq. is ap- 
proximately constant (see below). We used as least pa- 
rameters as possible to fit the data. The quality of the 
various fits is shown in the various figures that present 
the refractive index of He, H 2 , Ne, H 2 versus pressure. 



IV. DISCUSSION 

A. Polarizability 

The refractive index is generally related to the polar- 
izability of the molecular entities of the system through 
the the Lorentz-Lorenz formulation: 



LL = 



1 e- 1 



AnN A 



pe + 2 



(5) 



In Eq. Na is the Avogadro number, e the relative 
dielectric constant, and a the molecular polarizability 
of the compound. In the static case, this equation is 
called the Clausius-Mossoti relation. It arises from the 
expression of the local electrical field for dense materials: 
Eioc(r) = E 

macro 

(r) + 47rP(r)/3 + E'(r). P(r) is the 
polarization and E'(r) is the field created by the close 
neighbors of the electronic cloud centered on r. Eq. 
corresponds to E'(r) = 0. This is the case on symme- 
try center of a a cubic crystal. For hydrogen, helium 
and neon, Eq. 0is commonly usecUS^ because they are 
constituted by spherical molecular entities (in the case 
of hydrogen the molecules are freely rotating) , organized 
in highly symmetric structures. Also, for these systems 
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FIG. 6: Lorentz-Lorenz factor of H2O phases as a function of 
molar density. 
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FIG. 7: Lorentz-Lorenz functions of the studied compounds, 
normalized to their gas value-i^A 



in the pressure range covered here, the electronic transi- 
tions are far from the visible. The imaginary part of the 
refractive index should then be very small and e = n 2 is 
a very good approximation. 

In the case of water, the polarizability is generally 
written as the sum of three components: the polariz- 
ability due to the orientation of the permanent dipole 
of the water molecule, the polarizability due to the elec- 
tronic changes associated to the inter and intra- molec- 
ular modes and the polarizability associated to the elec- 
tronic transitions within the water molecules. The contri- 
bution of the permanent dipoles operates at frequencies 
below 10 5 Hz. The contribution due to the intermolecu- 
lar modes is important below 10 11 Hz in the liquid and 
10 5 Hz in the solid 22>. The contributions due to the 
intra-molecular modes should be important below 10 14 
Hz (the highest frequency vibration mode of H2O is a 
stretching mode at 3656 cm -1 in the gas phase). Above 
10 16 Hz electronic transitions should be considered. The 
optical window covered here is just in between these dis- 
persion phenomena, which implies a small dispersion of 
the refractive index and a small absorption. We did not 
observe any pressure effect on the dispersion. In Fig. 
the Lorentz-Lorenz factor of water and ice calculated us- 
ing the present data and published equations of state of 
liquid water, Ice VI and Ice VII (A JA, J&, M) is plot- 
ted, as a function of molar density p. It is interesting to 
note that the variation of LL plotted on Fig. El does not 
depend on the microscopic arrangement of the molecules 
but only of the density. This tends to prove that the 
hypothesis E'(r) = in the local field is also reasonably 
valid in water. Consequently, Eq. should be valid here 
to analyze the data. It is seen in figure 6 that the LL 
factor decreases by approximately 20 % up to the max- 
imum pressure of investigation. The approximation LL 
constant, that is used to analyse spectroscopic data of wa- 



ter under pressure, such as the Brillouin scattering sound 
velocity measurements—, is not correct because it leads 
to an underestimation of the refractive index of the order 
of 3 % at 35 GPa. It has been proposed from experimen- 
tal data below 1 GPa, that the evolution of the Lorentz- 
Lorenz factor of various molecular systems should follow 
an universal curve when plotted in reduced unitA, after 
small statistical dipole-induced dipole corrections. In our 
case, their relative effect on LL is smaller than 10~ 3 and 
can be neglected. In figured the change of the Lorentz- 
Lorenz factors of He, Ne, H2, H 2 vs. density is plot- 
ted in reduced units, i.e. as LL/LL gas (because LL gas , 
taken inSljlfiA contains the isolated atom polarizabil- 
ity) vs. Napct 3 , (where a is the Lennard-Jones parameter 
of the considered compound). The following Lennard- 
Jones parameters have been used: He, a — 2.556 
Ne, cr = 2.83 A^; H 2 , o = 3.06 A2i; H 2 0, a = 2.85 
Calculations have been performed using the follow- 
ing published EoS: HcAA Ne^A^A and H^AA- 
The scatter of the points give an idea of the experimen- 
tal errors. The large scatter at low pressure observed for 
He and H 2 , the most compressible materials, is mainly 
caused by the uncertainty on density calculation; more- 
over, uncertainty on LL is larger when n is close to one. 
At higher compression, relative uncertainty decreases. 

In Fig. a common trend is observed, namely a de- 
crease of the polarizability with density for the four sys- 
tems. In fact, the polarizability is related to the extent of 
the electronic cloud. In a dense environment, the atoms 
and molecules themselves compress to lower their inter- 
actions with their neighbors, and a decrease of the po- 
larizability is expected. Over the compression range cov- 
ered here, the scaled LL factors of H 2 and of He, two 
isoelectronic species, exhibit the same density evolution. 
However, in that pressure range, the fundamental differ- 
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FIG. 8: Refractive index vs. density of the studied com- 
pounds. The corresponded fitted linear functions are on ta- 
ble El 

ence between the density effect on electronic properties of 
these two media would be extracted from dispersion data. 
In He, it has been shown that the main contribution of 
the change of polarizability with density is taken into ac- 
count by the blue shift of the 1S-2P transition energy^. 
In contrast, the change of polarizability in hydrogen has 
been related to the red shift of the excitonic level X l Y,+ 
to B 1 ^^ that should lead to the metal hydrogenous,. 
But different evolutions are observed in Ne and H2O; the 
proposed universal behaviopiS. is nor valid nor theoreti- 
cally justified for a large perturbation by density. In the 
past 20 years, numerous computations and measurements 
have been performed for rare gas fluids (se o 41 ; 27 ) 42 and 
references therein), in order to reproduce, with success 
for rare-gas fluids, the two and three order refractivity 
virial coefficients. Yet, they have not been extended in 
the dense fluid or the solid phases. We are not aware of 
any calculation on the dielectric properties of water un- 
der pressure. On the other hand, the calculation of the 
dielectric properties of hydrogen is a subject of great cur- 
rent interest (see^ and references therein). We present 
here an interesting database of the optical dielectric prop- 
erties of low Z compounds under pressure. We hope that 
it could be helpful to validate ab-initio approaches. 

B. Gladstone-Dale relation 

The Gladstone-Dale relation expresses aa 4 ^: 



which leads to n = 1 + ap. It can be derived by differen- 
tiating Eq. [3and expanding it as a function of u = n— 1, 
if the polarizability is assumed to remain constant : 

du /„ 1 5 9 \ dp 
u \ o J p 



Compound 


a 


b (cm 3 /mol) 


He 


1 nnc _i_ n nm 
l.UUo ± U.UU1 


0. (49 ± 0.006 


r lulu IN e 


1 nnn -i- n nm 
l.UUU ± U.UU1 


1.410 ± u.uuy 


Solid Ne 


1.011 ±0.003 


1.25 ±0.02 


Fluid H 2 


0.994 ± 0.003 


3.26 ±0.03 


Solid H 2 


0.953 ± 0.007 


3.53 ±0.04 


Fluid H 2 


1.00 ±0.01 


6.05 ±0.2 


Ice VII 


1.081 ±0.01 


5.083 ±0.13 



TABLE II: Density linear evolution of the refractive index 
obtained in this work ; the index expresses as: n = a ± bp. 



Thus, this relation corresponds to the zero-order trun- 
cated series expansion of du vs. dp. For n < 1.8, the 
first order term remains smaller than 13 % of the zero- 
order term. Moreover, the volume dependency of a will 
tend to compensate the effect of the first and second or- 
der terms. 

The Gladstone-Dale relation has been widely used to 
analyze the shock- wave VISAR data. In fact, a trans- 
parent window is often attached to the rear face of the 
sample, together whith an optically reflecting interface, 
to prevent the formation of a reflected rarefaction wave. 
The high pressure refractive index of this window is 
needed to determine the rear face velocity of the sample 
by VISAR technique 4 ^. It can also been useful if the sam- 
ple itself is transparent, and if the shock front becomes 
reflective. This is the case for pre-compressed samples®. 
An advantage of the Gladstone-Dale relation is that re- 
fractive index expresses as a function of the density of the 
material alone (intrinsic temperature effect is not taken 
into account), and very simply. In this peculiar form, 
it allows to express straightforwardly the actual velocity 
of the reflecting interface as a function of its apparent 
velocity, for a planar single shock wave. 

For this reason, we think that a presentation of our re- 
sults in a form of Glastone-Dale fit could be useful for the 
high pressure community. They are presented in Fig. 03 
The data are close to straight lines, even if a small curva- 
ture is obvious in all curves. The Gladstone-Dale relation 
is often presented under the following form: n — a + bp^. 
We used this function to fit our data. The results are pre- 
sented on Table ITTI 



V. APPLICATIONS 
A. Equation of State 

The EoS of helium in the fluid and solid phase has been 
measured up to 13 GPa, as presented in Fig. EJ during 
runs 1 to 3. The volume of the sample was obtained by 
the multiplication of its thickness by its surface. Sample 
thickness was varying between 37 /im and 23 pva (simple 
gaskets). The evolution of the sample thickness during 
run 4 (double rhenium gasket) is also plotted in Fig. 
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FIG. 9: EoS of helium deduced from the sample thickness evo- 
lution upon pressure increase for runs 1, 2 and 3. It compares 
very well with previous determinations (— ™™- ), within 1 %, 
except for the lowest pressures. 
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FIG. 10: Evolution of the sample thickness upon pressure 
increase and decrease for Run 4. 



It is interesting to note that the thickness of the sam- 
ple decreased non-linearly upon pressure increase but re- 
mained approximately constant upon pressure decrease: 
from 47 fim to 48.5 /im. This can be used experimen- 
tally, when the sample thickness knowledge is needed^: 
one can perform measurements upon pressure decrease. 

Our EoS is compared with Brillouin spectroscopic mea- 
surements^ in the fluid phase and x-ray data in the solid 
phase^i^. At moderate pressures (P < 3 GPa), it seems 
that the compressibility of the sample is overestimated; 
this effect can be attributed to either leakage of the sam- 
ple or geometrical effects. Between 3 and 14 GPa, agree- 
ment with other data is excellent. No evident leakage 
has been observed in that pressure range. This shows 
that this method for EoS measurements is very inter- 
esting for liquids, where x-ray diffraction measurements 



are still difficult to perform and where there is no ra- 
dial pressure gradient in the sample. In particular, these 
measurements could be carried out for high temperature 
and high pressure liquids. 



B. Detection of phase changes 

The volume discontinuities at each phase changes can 
be deduced from measurement. In fact, we observed 
no evidence the the LL factor was discontinuous at any 
phase change studied. Thus, we assumed that this func- 
tion depends on density alone. By differentiating Eq. 
Ap/p can then be deduced from the observed An. The 
results are summarized tahle lllll Our data compare very 
well with those obtained by other techniques. Thus, 
careful refractive index measurement can detect phase 
changes with volume changes as small as 2 %. 

Phase change P (kbar) An Ap/p Ap/p 

xlO 3 (%) literature (%) 



Fusion He 


116 


3 


± 


2 


2 ± 1.3 




Fusion Ne 


47.7 


6 


± 


2 


3.5 ± 1.2 




Fusion H2 


53 


8 


± 


2 


3 ± 0.8 




Fusion H2O 


9.7 


41 


± 


4 


9.7 ± 1 


8.9^. 


Ice VI -> VII 


21.5 


37 


± 


4 


8.0 ± 0.8 


8.421 



TABLE III: Density discontinuities at phase changes; the tem- 
perature of measurements was between 297 and 299 K. 



VI. CONCLUSION 

We have performed accurate optical refractive index 
measurements of H 2 0, He, H 2 and Ne up to 35 GPa and 
at ambient temperature in a diamond anvil cell. These 
data should b useful to analyze some spectroscopic mea- 
surements on these systems at high pressure. The change 
of the refractive index has been expressed in terms of the 
change of the polarizability of the molecular entities in 
their dense surroundings. It is still challenging to com- 
pute such electronic properties. The present set of data 
could be useful in validating the various theoretical ap- 
proaches. Our method also allows the direct measure- 
ments of the Equations of State of transparent media: 
this has been carried out for fluid and solid He, with less 
than 1 % error. This method should then be interesting 
to determine the equation of state of dense fluids or dense 
amorphous media. 
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